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The C operator in PT-symmetric quantum mechanics satisfies a system of three simultaneous 
algebraic operator equations, C 2 = 1, [C,PT] = 0, and [C,H] = 0. These equations are difficult 
to solve exactly, so perturbative methods have been used in the past to calculate C. The usual 
approach has been to express the Hamiltonian as H = Ho + eHi, and to seek a solution for C in 
the form C = e®V, where Q = Q(q,p) is odd in the momentum p, even in the coordinate q, and 
has a perturbation expansion of the form Q = eQi + e 3 Qs + t'Qs + . . .. [In previous work it has 
always been assumed that the coefficients of even powers of e in this expansion would be absent 
because their presence would violate the condition that Q(p, q) is odd in p.] In an earlier paper it 
was argued that the C operator is not unique because the perturbation coefficient Qi is nonunique. 
Here, the nonuniqueness of C is demonstrated at a more fundamental level: It is shown that the 
perturbation expansion for Q actually has the more general form Q — Qo + eQi + e 2 Q,2 + ■ ■ ■ in which 
all powers and not just odd powers of e appear. For the case in which Ho is the harmonic-oscillator 
Hamiltonian, Qo is calculated exactly and in closed form and it is shown explicitly to be nonunique. 
The results are verified by using powerful summation procedures based on analytic continuation. It 
is also shown how to calculate the higher coefficients in the perturbation series for Q. 

PACS numbers: 11.30.Er, 03.65.Fd, 02.30.Mv, ll.10.Lm 



I. INTRODUCTION 



The properties of PT-symmetric Hamiltonians have been observed in a wide variety of laboratory experiments 
[Tl410j . For a PT-symmetric Hamiltonian having an unbroken PT symmetry a linear 'PT-symmetric operator C exists 
that obeys the following three algebraic equations: 

C 2 = 1, (1) 
[C,VT\ = 0, (2) 
[C,H] = 0. (3) 

Constructing the C operator is the key step in showing that time evolution for a non-Hermitian PT-symmetric 
Hamiltonian is unitary [TTl [T2"] . 

The C operator for a few nontrivial quantum-mechanical models has been calculated exactly |13H17j by solving 
Q-Q. However, in general this system of equations is extremely difficult to solve analytically. Therefore, in most 
cases a perturbative approach has been adopted for the solution of these equations fT8H2"2"] . 

The standard approach to solving ([l])-(|3]) has been to express the C operator in a simple and natural form as an 
exponential of a Dirac Hermitian operator Q multiplying the parity operator P: 

C = e Q V. (4) 

Note that e^/ 2 is precisely the metric operator r\ discussed in Refs. [2"3Tf2"8] . This operator can be used to construct a 
similarity transformation that maps the non-Hermitian Hamiltonian H to an isospectral Hermitian Hamiltonian 29 . 
If we seek a solution for C in the form Q , we find that (JlJ and Q , which can be thought of as kinematical equations 
because they hold for all choices of H, imply that Q(p, q) is an odd function of the momentum operator p and an 
even function of the coordinate operator q [TTJ Q2] . The problem is then reduced to finding the solution to ^ , which 
can be thought of as a dynamical equation because it refers to the Hamiltonian H . 
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It is difficult to find a closed- form analytical solution to ([3]). However, in the past this equation has been solved 
perturbatively as follows: Express the Hamiltonian in the form H = Hq + tH\ and treat e as a small parameter. 
Then, seek Q as a formal perturbation series in odd powers of e: 

oo 

Q(p, (? ) = ^e 2 ^ 1 Q 2j+1 (p, g ). (5) 

3=0 

The obvious question to ask is, Why do only odd powers in e appear in the perturbation series The explanation 
that has been given in the past is that even powers of e are excluded from the series because Q(p, q) is required to be 
odd in p and even in q. The reasoning goes as follows: In the quantum-mechanical cases that have been studied so far, 
such as H = Hq + ieq and H = Hq + ieq 3 , the unperturbed Hamiltonian Hq = + \q 2 is the harmonic-oscillator 
Hamiltonian. If there were a term Qq€° in the series jHJ, then Qq would satisfy the commutation relation [Qq, Hq] = 0. 
The vanishing of this commutator implies that Qq is a function of Hq, and thus it is an even function of p, which 
shows that Q = 0. Once it is established that Qq = 0, it is relatively easy to show (see Ref. [TS], for example) that 
Qij — (J = 1, 2, 3, . . .). We show in this paper that this argument is actually incorrect; there are indeed solutions 
to the commutator equation [Qq,Hq] = that are odd in p - infinitely many such solutions, in fact. It is precisely 
because of the existence of these odd-p solutions that the C operator is nonunique. 

It has recently become clear that the nonuniqueness of the C operator has important implications for the mathe- 
matical and physical interpretation of 'P7'-symmetric quantum mechanics [30 - 32J . In Ref. [31] it is shown that if the 
C operator is nonunique, then it is unbounded, and in this paper we verify this result explicitly. 

The paper |33j is relevant because it discusses for the first time the existence of multiple (nonunique) solutions to the 
commutator equation ([3]). In Ref. |33| it is shown that the inhomogeneous commutator equation [Q\,Hq] = 2Hi has 
an infinite number of particular solutions, which differ from one another by solutions to the associated homogeneous 
commutator equation [X, Hq] — 0, and it recognizes that the solutions X are not necessarily functions of Ho only. 
The importance of nonunique solutions to commutator equations is also central to Ref. |34j , where a particular time- 
operator solution O to the inhomogeneous commutator equation [Q, H] = i is called minimal and a classification of 
the infinite number of associated nonminimal solutions is given. 

The approach used in the current paper is based on the recognition that finding multiple solutions for Q\ in ([5]) is 
not the only way to demonstrate nonuniqueness. Here, we introduce a clearer and more fundamental way to explain 
the nonuniqueness of the C operator. We show that a more general way to represent Q is by the expansion 

oo 

Q{p,q)=J2^Qj(p,<l) (6) 
j=o 

in which all nonnegative integer powers of e appear. An advantage of this new representation is that in the limit e — > 
we obtain an infinite class of exact C operators for the quantum harmonic-oscillator Hamiltonian Hq = ^p 2 + \q 2 ■ 
Then, once we have Qq for the harmonic-oscillator case, we can straightforwardly generalize this result and verify 
that the operator C is nonunique. 

This paper is organized as follows: In Sec. [H] wc show how to construct exact and explicit closed-form solutions 
that are odd in p and even in q to the homogeneous commutator equation [Qq(p, q), Hq] = 0. Our result is that there 
is a unique bounded C operator and a nonunique infinite class of unbounded C operators for the quantum harmonic 
oscillator. The techniques used in Sec. [TT] involve the formal summation of infinite series of singular operators. However, 



in Sec. Ill we verify the validity of the formal calculations done in Sec. [H] by applying powerful summation techniques 
that are used to regulate divergent Feynman integrals. This verification leads us to conjecture that it may be possible 
to apply the principles of summation theory to extend and generalize the rigorous notions of Cauchy sequences and 
completeness expansions, which are used in mathematical Hilbert-space theory, to divergent sequences and series of 



vectors. Next, in Sec. IV we develop the formal machinery needed to determine the higher coefficients Qi, Q2, ■ ■ ., 
in the expansion ([6]), and in Sec. |V]we concentrate on calculating Qi(p,q) for the specific case Hi — iq. Finally, in 
Sec. |VI| we make some brief concluding remarks. 



II. SOLUTIONS TO [Q (p, q), H ] = THAT ARE ODD IN p AND EVEN IN q 

A powerful strategy for solving operator equations of the form 

[Q o (p,q),H o ] = (7) 

is to represent the solution Qo(p,q) as an infinite series of totally symmetric operator basis functions T m n . The 
operators T m>rl are described in detail in Refs. |34H37j . However, to make the presentation in this paper self-contained, 
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we recall that for m, n > the operator T m n is defined as a symmetric average over all orderings of m factors of p 
and n factors of q: 

T i,i = l(pq + qp), 

T\,2 = l(pqq + qpq + qqp), 
and so on. The operators T m ^ n obey simple commutation and anticommutation relations: 





ifiT rnn _ \ , 
















= — 2mT m _(_x 5 n-i 




= 2zmT rn _i jn _(_i, 



(8) 

where the curly brackets indicate anticommutators. The operator T mj „ can be re-expressed in Weyl-ordered form 36J: 



A:=0 V 7 fc=0 V 7 



p m q n ~\ (9) 



where m, ti = 0, 1, 2, 3, ■ • • . Introducing the Weyl-ordered form of T m n allows one to extend the operators T m n 
either to negative values of n by using the first sum or to negative values of m by using the second sum. The 
commutation and anticommutation relations in Q remain valid when m is negative or when n is negative. 

To find solutions that are odd in p and even in q to the commutator equation ([7]), we take Qo(p,q) to have the 
general form 

^ 7 W = E a i 7) ^ 7+ i-2 fc ,2 fc , (10) 

k 

where 7 = 0, ±1, ±2, ... is a parameter. Substituting (10) into (|7|), we obtain the following two-term recursion 
relation for the coefficients : 

4 7 +i(fc + l)-(7-fc + l/2)4 7) -0 (k = 0,1,2,...)- (11) 

This recursion relation is self-terminating; that is, if we choose a_l = 0, then is an arbitrary constant, aj, 7 " 1 
vanishes for k < 0, and a^ 7 ' for k > is determined in terms of as the solution to the recursion relation 



(* = 0,1,V..). ,12, 



The series (10) with coefficients (12) can be summed as a binomial expansion: 

oo 

^a^^a^l + x 2 )^ 2 (13) 



Thus, for each 7 > the odd-p and even-g one-parameter family of solutions to the homogeneous commutator equation 
is 

< * , -^["{{( i+ «Hr' 2+ ( i+ H«r''-44"-^ • <i4) 



where we have used the identity [33 

T „„=,.. , 

2 V pj 2 \q 



(27+I) times 



T-n,n = Uq-) +1(-P) ■ <>">) 
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As stated in Sec.|lJ we can see that while commutes with the Hamiltonian Hq — \p 2 + \q 2 , it is not a function 
of Hq because by construction it is odd in p. Furthermore, while the construction of the solutions in ( 14 ) involves 
scries in inverse powers of p, these solutions are well behaved as p — > 0. To sec explicitly the oddness in p we display 
the solution corresponding to 7 = 0: 



,(0)1(0) I L 1 l /..ll /.il /.il 



Qo = 7«0 \ l + q-q- P+P\ l + q-q- + \ 1 + -q-q P + P \ 1 + -q~q ■ (16) 
4\V P P V P P V P P V P P J 

In the classical limit for which p and q become commuting numbers, this solution becomes 

Classical ^SgnWV^T^; (17) 

the oddness in p is evident. 

It is important to point out that for the harmonic oscillator, which corresponds to e — 0, the metric operator 

77 = is just unity when Qo — 0. Thus, the metric operator is bounded of this special case. However, for Qq in 



(16 1, the metric operator is no longer bounded, but rather for large q it behaves like e q and for large p it behaves like 
e p . Needless to say, since there is an infinite number of possible choices for Qo, there is an infinite number of possible 
metric operators. Only one of the metric operators is bounded 



III. USING SUMMATION TECHNIQUES TO VERIFY RESULTS OF SEC. [XT] 



We observed in SecOthat even though solutions Qo 



(7) 



in ( 14 ) were constructed by performing a formal infinite sum 
over arbitrary powers of the inverse momentum operator 1/p, these solutions are well behaved as p — ¥ 0. However, the 
calculations in Sec. [H] are certainly not rigorous. The aim of this section is to provide mathematical support for the 
validity of the formulas in (14 1. Specifically, since the C operator commutes with the Hamiltonian, the nth eigenstate 
\ip n ) of the Hamiltonian must also be an eigenstate of C. We expect that the eigenvalue of \ip n ) is (—1)™. For this to 
be true, \ip n ) must be an eigenstate of Q ^ with eigenvalue for all n: 



? 7) hM = o. 



(18) 



In this section we show by explicit calculation in which we use powerful summation techniques that this is indeed the 
case. Here, we limit our calculation to the case 7 = 0. We first consider the ground state \i)}q) and then generalize to 
the nth eigenstate. 

From (10) and (12) we see that Q °^ is given by 



}(0) 



= ^2 a fc Ti_ 



2fe,2fc) 



a k = a - 



k=0 



(-l) k T(k- 1/2) 



(19) 



Also, the unnormalized eigenfunctions of th e harmonic-oscilla tor Hamiltonian in coordinate space are given by ipn(q) 
H n {q)e~ q l\ While the formal sum in il9l can be written as (U6|, it is difficult to use this result to verify the eigenvalue 



equation ( 18 



to perform t 



. A better strategy is to calculate the action of each term in the sum in ( 19 ) on the eigenstates and then 
le summation over k. 



Because inverse powers of the momentum operator arise in (19), it is most convenient to work in the momentum 



representation, where the eigenvalue equation ( 18 ) becomes 



(pIQ 0) I^) 



1 00 



a fc (-l) A 



fe=0 



p { -d, 



2k- 



P 9 ; 



2/. 



21: 



2k 



P \ ^nip) = 0, 



(20) 



where we have used (15) and the third formula in 
of t/j. 



81). Here, ip n (p) = (-j)"^n(p), where ip is the Fourier transform 
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A. Ground state tpo 



Let we first study the simplest case n = 0. We can see that the action of each term in the series expansion (201 on 
the eigenstate tpo produces more and more negative powers of the momentum p: 

n( o) M r M If, /, 2\ / 12 48 90 \ / 30 240 1350 5040 \ 

Q { ) (P)MP) = -|4p ao -^4 P +-ja 1+ ^+- + - + ^Ja 2 -^+- + — + — + — j a 3 



56 672 6300 47040 264600 997920 1891890 \ , „*, 2 

4p H — -\ = — | =- -\ — H 1 1 — a 4 - . . . > e p /2 

pO pO pi pi> pii pia pxo j [ 



We can rearrange the series in (21 1 to read 



r oo oo r_-.sk 

I fc=0 fc=0 ^ 



where P n {y) are polynomials: 



(21) 



Q<%)Mp) = | <| 4p^(-l) fe a fc - 5] i_J_ afc+1 P 2fc(2| ,2 ) I, e -P 2 /2 ; (22) 



^(,) = E 2 iV^W +2) V - <»> 
2 Q a! a + 2)! n — a ! 

To verify that Qq (p)V>o(p) = 0, we must show that 

^ r(fc-l/2) 1 A 1 (2a + 2)! ~ r(2fc + l)r(fc + 3/2) 
P ^ fc! p 3 ^ p 2Q a!(a + 2)!2« ^ T(fc + l)r(2fc - a + 1) ' 1 J 

where we have substituted the above formulas for the polynomials P n and the coefficients a^. It is easy to verify 
that the exact sum of the convergent series XlfeLo-^(^ — l/2)/fc! is zero. However, the series X^fcLo r(k+tyrf2k-a+i) * s 
divergent for a > 3/2. Therefore, we must introduce a summation procedure to make sense of this series. 

Our summation procedure is a discrete variant of dimensional continuation, a technique that is used to interpret 
divergent Feynman integrals. To illustrate our approach, let us consider the following D-dimensional integral: 

m=ld D x^±^ 2 . (25) 




FIG. 1: A plot of g(D) in p6l for < D < 6. Note that g(3) vanishes even though the integrand of 1(D) in (25| is strictly 
positive for all D. This shows that the notions of positivity and negativity evaporate in the case of a divergent integral 
representation. 



() 



This integral converges for D < 2 and its exact value is 1(D) = Sd g(D), where Sjj — 2it D / 2 /T(D/2) is the surface 
area of a D-dimensional sphere of radius 1 and 

Evidently, even though the integral representation for 1(D) diverges for D > 2, the function 1(D) is well defined and 
analytic in D except for isolated simple poles. [Figure [I] gives a plot of g(D) for < D < 6.] Observe that 1(D) 
vanishes for D = 3. This result is surprising and somewhat counterintuitive because the integrand of 1(D) is strictly 
positive when D — 3. The vanishing of 7(3) shows that when an integral representation is divergent we cannot draw 
qualitative conclusions regarding the sign of its value. (Indeed, the Borel sum of the series 1 + 2 + 4 + 8 + .. . is 
uniquely —1 even though all of the terms in this divergent series are positive!) 



We will now show that while the sum over k in (24 1 diverges for a > —3/2, we can evaluate the sum for a < —3/2 
and then use analytic continuation in a to sum the series for a = 0, 1, 2, 3, . . .. The surprising and counterintuitive 
result is that while the summand is positive for 2k+l > a, the sum of the (divergent) series vanishes for all nonnegative 
integer values of a. 

The divergent series to be summed is 

_^ r(2fc + i)r(fc + 3/2) 

" 2^T(k + l)T(2k-a + iy [ ' 



To perform the sum we first express the coefficients in terms of the beta function: 

r(a;)r(y) 
T(x + y) ' 



= 3^. (28) 



whose integral representation is 

B(x,y)= f dtt^' 1 ^ (l-tfy-V (Re(ar), Re(y) > 0). (29) 
Jo 

By multiplying and dividing F a by T(—a) we obtain the result 

1 v£(fc + 3/2) 1 - T(fc + 3/2) f* 2k t 

Fa -f(=a)2^ T(k + 1) B ^ k + l ^-T^a)^ n T(k + 1) //" (1 " <} " (30) 



k=0 v ' v ' fc=0 

We then use the binomial expansion 



^ r(fc + 3/2) 2fc _ VtF _ ,2w3/2 

I-* k\ - 2 [ ' 



k=0 



(31) 



to show that 



^ L 1 - (r ^ 5/2(1 + tr " 2 - ^ r irT- a ") 3/2) - Fi {>■- - H ■ (32) 

This function vanishes for all nonnegative integer values of a because the hypergeometric function is finite for these 
values of a. In Fig. [2] we plot F a for —5 < a < 10. Observe that F a vanishes for a = 0, 1, 2, Note also that F a 



is singular at a = —3/2, the value of a for which the series (27 1 begins to diverge. Interestingly, this function has 
double poles at the half-odd integers a = 1/2, 3/2, 5/2, .... As a increases, the double poles begin to resemble single 
poles 



B. Generalization: eigenvalue equation for ip n 

In this subsection we study the action of the operator on the nth eigenstates of the harmonic-oscillator 

Hamiltonian. We describe first the case for which n is even. (The odd-n case is treated in an identical fashion and is 
considered briefly at the end of this subsection.) For even n 



OO / 

E 

k=i 1 



-l) fe , 



E 

3=0 



p 23 + iST(-l) k a k S^(k) 



k=0 



V/2 



(n= 1,2,3,...). 



(33) 
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FIG. 2: A plot of F a in Q for -5 < a < 10. 
r(— a) in the denominator is infinite when a — 0, 



Note that F a vanishes at every nonnegative integer. This happens because 
1,2,... and the hypergeometric function is finite. These zeros are all simple 



zeros, but the poles at the half-odd integers beginning with 1/2 are all double poles, as is verified in Fig. [3] 




FIG. 3: A blow-up of the region near a = 9.5 in Fig. [2] In the left panel is a plot of (a — 9.5)F a for 9.496 < a < 9.505 and in 
the right panel is a plot of (a - 9.5) 2 F a for the same range of a. Note that the left plot has a typical simple-pole behavior at 
a = 9.5 while the graph in the right panel is finite at a = 9.5. 



Both the sums over k in (33) diverge for a = 0, 1,2, . . . except for j = n in the first series. (The special case j = n 



gives the convergent series that was already considered in Subsec. Ill A ) The divergent series can be evaluated by 
using the summation procedure that was introduced in the previous subsection. Following the procedure for the n = 
case, we will show that if we write ([33]) as 



Q^2n 



3=0 



,-p 2 /2 



(34) 



then both M (n ^ and A/" (n) vanish for n = 1, 2, . . ., j = 0, 1, . . .n. 
We first consider the scries 



oo / -i \k 



(35) 



k=l 



where j£ (x) are fcth order polynomials that can be written as 



v 7r (« — a)!(n + a + 1)! 

and E n ^ a (k) are polynomials of degree n in the variable k. The first four polynomials are 



(36) 



10a, 



E lj0 (fc) = 2(1 - a)fc + a 2 - 4a, 

E 2>a (k) = (4a-8)fc 2 -(4a 2 -20a + 4)fc + a 3 -7a 2 
£ 3 ,a(fc) = (24 - 8a)fc 3 + (12a 2 - 72a + 24)fc 2 - (6a 2 - 48a 2 + 70a - 24)k + a 4 - 9a 3 + 20a 2 - 48a, 
E i>ol {k) = (16a- 64)fc 4 - (32a 2 - 224a + 96)/fc 3 + (24a 3 - 216a 2 + 320a- 224)fc 2 

-(8a 4 - 80a 3 + 176a 2 - 536a + 96)fc + a 5 - 10a 4 + 23a 3 - 158a 2 + 216a. 



In terms of E n>a (2k — 1) = X) 7 =i e a, 7 fc 7 7 the first term in (34) can be written as 



(37) 



A^) = Iran - D" T 2 " r ( a + 1 / 2 ) y 

yV "tt 1 J " 2^ „2a+l r« . + rv + 1)! 2^ 



where A/" a , 7 is given by 



n p"*-<- x (n + a + 1). 

Q— 7—1 



_ - r(fc-i)r(2fc + i) 



A/" a , 7 , 



(38) 



(39) 



This is the divergent series that we need to study. 

Because 7V Q , 7 is divergent for a > —5/2 — 7, we evaluate the sum for a < —5/2 — 7 and use analytic continuation in a 
to sum the series for a = 0, 1, 2, . . .. Multiplying and dividing jVq., 7 by T(— a — 1) and using the integral representation 
of the beta function B, we obtain 



T(-a- 



1 „1 oo 



The sum over k in ( 40 ) gives 



fe=0 



r(fc) 



fc=0 



7+1^7 ( 2' 2 ' 



r(fc + i/2)fc^ 2fc 
r(fc) 



(40) 



(41) 



where the first dots in the hypergeometric functions stand for (7— l)-twos and the other dots stand for (7 — l)-ones. For 
fixed 7, the hypergeometric function in (41 ) can be written as I/ 7 (i 2 )2~ 7 (l — t 2 )~ 3 / 2 ~ 7 , where £ 7 (i 2 ) = Y11=i ^, 7 i 2<T 
is a polynomial of order 7 in the variable t 2 . In terms of L 1 {t 2 ) the series (39) becomes jV Q)7 = Y^Z=i ^<r,yJ^a,-Y,<r: 
where 



1 



7,(T 



2 T(-a-l 
% /7rS(5 + 2cr,-a-7-f 
~2 r(-a - 1) 



dtt 4+2CT (l - i)- Q - 7 - 7 / 2 (l + iJ-i-T 



3 5 
-2-F1 7+ -,5 + 2cr;-a-7 + 2o- + -;-l 



(42) 



The function J\f a ,7,<r vanishes for all nonnegative integers a, 7, and a because the denominator becomes infinite while 
the hypergeometric function is finite. The special case M a .\,\ is plotted as a function of a in Figs. [4] and Fig. [5] The 
vanishing of M a ^,a guarantees that the first sum in (33) is identically zero. 

Next, we evaluate the sum over k in the divergent series in the second term in (33): 

00 

M = £(-i)Wft(fc), 



(43) 



k=0 



9 



ft 



-6 



-4 



-2 



-4 



-2 



FIG. 4: A plot of jV a ,i,i in (421 for —6 < a < 8. Like the special case displayed in Fig. [2] Af a ,i,i vanishes for all nonnegative 
integers. However, A/" a ,i,i is different from F a in that it has simple poles rather than double poles. It is not completely obvious 
that A/" Q ,i,i vanishes when a = 0, 1, 2, . . ., so in Fig. [5]the behavior of 7V Q ,i,i near a — 5 is blown up. 



0.15 



0.10 



0.05 




4.95 



"5.00 



5.05 



5.10 



FIG. 5: A blow up of the graph of Af a ,i.i in Fig. [4] near a = 5. Observe that there is a zero at exactly a = 5 as well as an 
additional zero near and to the left of a — 5. 



where S^.(fc) = E£ -o Snjk 1 are polynomials of degree n — j in the variable k. The polynomials are listed below for 
n = 1,2,3: 



S^ik) = 8(8fc + l), S^\k) = W, 



S?\k) = 16(68fc 2 + 14fc + 3), S\*'(k) = -64(8fe + 3), S*^(fc) = 64, 



(2), 



Sf^/c) = 96(192fc 2 + 52fc 2 + 46£; + 5), Sg J (k) = -64(196fc 2 + 142fc + 45), 



(3)/ 



Sfp ) (Jb) = 384(8A + 5) 1 ST( fc ) 



(3), 



-256. 



(44) 



As noted above, the series for j = n [the highest power in p in (33)] are convergent and their exact sum is zero. For 



the divergent cases, the series (43) can be written as M — J2l=o s n,iMt, where Mi is the divergent series 



k=0 



fc! 



(45) 
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This series can be rewritten in the form 



V ^ Jo i — n 



(46) 



fc=0 

Summing the series in (46) over k, we obtain Y^k=Q^ k ^ = (1 — X^=o r t-<?t a ■ Thus, A4( can be written as 
•M-l = Y,i=o r ^^M^cr, where 



M a e = - 



1 



dt(l - t )-5/2-^-3/2+ CT 



2^ 



B(<7- 1/2,-^-3/2), 



which vanishes at zero for (a, £) = 0, 1,2... .. 

Finally, we consider the case of odd n. For this case we get 



2n+l 



n+1 oo 



^4feT^ a fe°2fc (P ) 



j=0 k=0 



k=0 



V/2 



Here, V^["^(fc) are polynomials of degree n — j in the variable k: 
Vi (1) (fc) = 16*, V (1) (/c) = -8, 



v¥\k) = 16(14fc 2 -fc), F^Cfe) = -48(4fc + l), V^>(k) = 32, 



,(2) 



(2), 



V 3 (3) (fe) = 128(26fc 3 -4fc 2 + 3fc), V^'(k) = -32(124fc 2 + 58/c + 15) 
^ (3) (fc) = 640(2/c + 1), V (3 \k) = -128. 



(3), 



Also, the polynomials (x) can be written as 

2™+ 3 



0Wte) - 2 —(2n + 1)" V ^(a + 3/2)(fc + 2)! (a) 
^ fc W v^tF 1 (fc-a)!(n + a + 2)! ^ n l j 



where Un (k) are polynomials of degree n in the variable fc. The first four such polynomials are 

U[ a \k) = 1, 

U ( 2 a \k) = 2/c-a + 3, 

C/^ a) (fc) = Ak 2 - (4a - 12)fc + a 2 - 56a + 12, 

Ui a) (k) = 8/c 3 - (12a - 36)/c 2 + (6a 2 - 30a + 76)fc - a 3 + 6a 2 - 29a + 60. 



(47) 



(48) 



(49) 



(50) 



(51) 



As in the case of even n, we can verify that the sum of both series in (48) vanish for all n. 

This completes the verification that the eigenfunction i\) n of the harmonic-oscillator Hamiltonian is also an eigen- 
function of the Qq operator with eigenvalue 0. Thus, ip n is an eigenfunction of C = eQ°V with eigenvalue (—!)"■ 



IV. CALCULATION OF Q TO FIRST ORDER IN e 

The general approach in this paper is to find an operator C = e®V for the PT-symmetric Hamiltonian H = Hq+cHi, 
where Q has the power-series expansion ^ in the parameter e and this expansion has a nonvanishing zeroth-order 
term; that is, Qq ^ 0. In this section we concentrate on the formal problem of determining the first-order coefficient 
Qi once Q is given. 

The coefficient Q\ satisfies the equation 

[e Qo+eQl ,H ] =e{e Qo ,Hi}, (52) 
which follows immediately from ([3]). If we expand this equation to first order in e, we obtain 

Z + l -(Q a Z + ZQ ) + l -(Q 2 Z + Q ZQo + ZQ 2 ) + ^(Q 3 Z + Q 2 ZQ Q + Q Q ZQ 2 + ZQ\) + ...= {e % , i^} , (53) 
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where 

Z=[Qi,ff ]. (54) 

Recall that Qq is a solution to [Qq,Hq] = 0, which is a homogeneous equation. Thus, any parameter /x times Qo 
is also a solution. Our approach will now be to make the replacement Qo ~~ ► fJ-Qo i n 1 53 ) and to treat \i as a small 
perturbation parameter. We can thus expand Z as 



Z = J2 Z ^ n (55) 



To zeroth order in /i we obtain the result Z$ = 2Hi. To first order in fi we obtain Z\ = 0, and in fact we find that 
Z^j+i = for j = 0, 1, 2, . . .. The general result for n > 2 can be given in terms of Bernoulli numbers B n : 

2B 

Z n = [Qo, ■ • • [Qo, [Qo,Hx]} . . .]„ timos , (56) 
where Bo = 1 > B\ — ~ \ (which is not used in the above formula) , and 

B 2 = \,Bz = 0, Bi = B 5 = 0, B 6 = £ 7 = 0, £ 8 = -i, £ 9 = 0, B w = A, Bn = 0, B 12 (,< " 



6' J ' " 30' u ' u 42' ' ' ° 30' » ' ±u 66' ^ ' " 2730 
We now decompose Qi into a perturbation series in powers of ji, 



DC 



Qi = ]T R ^ n > ( 5? ) 



7i = 



and obtain a sequence of commutator equations for the coefficients R n : 

[R n ,H Q ] = Z n (n = 0,1,2,...). (58) 
In the next section we show how to solve these commutator equations for the special simple case in which Hi = q. 



V. SOLUTION OF d58j) FOR THE SHIFTED HARMONIC OSCILLATOR H = \p 2 + \q 2 + ieq 



Let us consider the shifted harmonic oscillator for which Hq = \p 2 + \q 2 and H\ = iq. This Hamiltonian has an 
unbroken VT symmetry for all real e. Its eigenvalues E n = n + \ + \t 1 (n = 0, 1, 2, . . .) are all real. One C operator 
for this theory is given exactly by [11] [12] 



C = e~ 2ep P. 



(59) 



In the limit e — > the Hamiltonian becomes Hermitian and C in (59) becomes identical with V. However, the solution 
for C in (59) is not unique, and by taking any or all of the Qo in (14), we obtain an infinite number of operators C. 
To find Q\ we must calculate Zq (which is 2iq), Z 2 , Z4, and so on, and from these we must solve (58) to obtain i?o, 
i?2j Ri, and so on. 



For the case n = in (58) we have a simple exact solution to the commutator equation for Rq: 

Ro = -2p. 



(60) 



We emphasize that this solution is not unique. 
The equation for Z 2 is 



00 00 



4EE4 7) «i 7) d-2fc-2 7 )4i\ 



3=0 k=l 



(61) 



where F, 



(7) 



[T_ 



■■j t k - L J --2j-2 7 +l,2i,7 1 _2fc-2 7 

elements T m n : 



2k], whose explicit form is obtained by using the algebra in Ref. [33j of the basis 



j+k-l 2a+l 

F ji = E E CjX,a,/3 :r - 2 ( fe +j)- 4 7-2a,2(fc+i)-2Q-l 
a=0 0=0 



(62) 
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with coefficients C a g given by 



C 



(7) 



i a (-l) a+ P(2j)\(2k)\(2j + (3- 2)!(2fc + 2a - 0)\ 



j,k,a,P 4«(2j - 2)!(2fc - l)!(2j + $ - 2a - l)!(2fc - /3)\(2a + 1 - /?)!£! ' 



A more compact form for Z2 in (61 1 is 



00 00 



Z 2 - -g X! X! ^fcTa T -2fc-4 7 -4a,2fc-l, 
fc=l a=0 



where 



(7) = (-l) fc -« +1 r(fc + 2a + 2 7 )r 2 (a + 7 + 1/2) 

k < a r( 7 - l/2)r(fe)r(a + i)r(a + 27) (a + 7) 



(63) 



(64) 



(65) 



and for simplicity we have set ao = 1. 

For the special case 7 = we have the following results for i? 2 °\ where from now on we omit the superscript (0). 
The commutator equation is 



oo oo 



(66) 



k=l a=0 



This is a linear equation, so we solve it for each a separately and express the solution as a sum over a: i?2 
Sa=o Pk,aR2,a- For a = we seek a solution of the form 



R2.0 — ^ PkflT-2k-l 



2k 



fc=0 



whose coefficients p/^o satisfy the recursion relation 

(2k - l)/3fc-i,o + 2kp k ,o = A kfi . 



(67) 



(68) 



(The techniques used here are described in detail in Ref. [33] .) The simplest solution to this recursion relation is 

Pkfi = 7T(-l) fe . 

For a = 1 we set 



R2.1 — ^ Pk,\T-2k-h 



■11; 



fe=0 



and so the recursion relation for the coefficients p k) i is 

(2k + 3)p k -i,i + 2kp kfi = A kA , 

whose solution is p k ,o — — f (— l) k (k + 2)\/k\. 
For general a we have 

oo 

R2,a — Pk,gT-2k-ia-l,2k 
k=0 

and the recursion relation for the coefficients p k , a is 

(2k + 4a - l)pk-i, a + 2kp k:Ct = A k , a , 

whose solution is 



(69) 



(70) 



Pk,, 



k (k + 2al 
K ' k\ 



r(a + l/2) 



(71) 



(72) 



(73) 
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A. Complete evaluation of the first-order expansion Q\ 

We now derive the general form of the first-order expansion in e of the operator Q = Qo + eQi + e 2 Q 2 + ■ ■ ., which 
takes the form of a series in even powers of the parameter \x: 

oo 

Q 1 = J2» 2nR 2n, (74) 

n=0 

whose coefficients R 2n satisfy the commutator equation 

[R 2n ,H ] = Z 2n (n= 1,2,3,...). (75) 

Recall that the operator Z 2n is proportional to the recursive evaluation of the double commutator [Qo, [Qo, • ]] 
acting on the operator Z 2n - 2 , starting from Z = 2H\. (As established earlier, Z 2 — [Qo,[Qo,Hx]]/6, Z4 = 
~ [Qo, [Qo, [Qo, [Qo, -Hi]]]]/360, and so on.) The operator Z 2n can be written as a double series over the basis el- 
ements T miTl once their algebra has been repeatedly applied and its closed-form expression is 

00 00 

Z 2n = ^ A< k2 )T -2k-4a,2k-2n+l (n = 1,2,3, .. .), (76) 

k—n a—0 

where the coefficients A^ 1 ^ are given by 



= Wi 2n \-l) k T(k + 2a + 7i-l)/T(k). (77) 

ous to evalu 

for n — 1 we have Wa^ = [T(a + l/2)/a!] 2 ; the evaluation of the next commutator with Q gives 



It is extremely laborious to evaluate explicitly the function W a 2n \ even for the first few values of n. For example, 



(2) = " T(2a -(3 + lJTjp + l/2)r 2 (a -0+ 1/2) 
h'o ( 4a + ^PW&a -(3 + l/2)r 2 (a - /? + 1) 



Moreover, the existence of solutions R 2n to ( 75 ) is not affected by the explicit form of the functions W a 2n ^ because 
(75) is a linear equation. 



Substituting (76) into (75) and noting the last two commutator equations in (|8|), we argue that the operator R 2n 
has the form 



oo oo 



R%n — Pfc iQ ' > 7 1 -2fe-4Q-l,2fc-2n+2, (78) 



k=n— 1 a—0 



where the coefficients p^ 2 ™'* satisfy the recursion relation 



{2k + ±a- l)p£3 ia + 2(k - n + l)p£? = A<j£ . (79) 



Equation ( 79 ) is a first-order linear inhomogeneous difference equation that we can rewrite as 

(2 „) 2k + 4q + 1 (n) _ ^k+la rsn s 
P *+i><* + 2(k-n + 2) Pk > a 2(fc-n + 2)' 1 ' 



To find solutions to ( 80 1 we divide both sides of the equation by the summing factor Yk , where 

k 



y = ( i \ k TT (2j ± Aa ± 1} = r(fc + 2a + 3/2) r>m 

k 1 J 11 2(7- n + 2) r(2a + n + l/2)r(/c-n + 3)' 1 ' 



Equation ( 80 ) then becomes 

Pk + l,a Pit, a A k+l,a 



„(2n) J2n) A (2n) 



Y k Ffe_i 2(fc - n + 2)Y k ' 



(82) 
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Note that (82 1 has taken the form of an exact discrete difference of the function pi 2 ™ /-Pfc-i- Summing both sides of 

G (2n) 



( 82 1 from 1 to k — 1 gives the solution 

(-l) fc+1 r(fc + 2a + l/2) 
r(n + 2a + l/2)(fc-ri + l) 



(2n) 
Pk,a 



k-1 



3=1 



r(n + 2a + l/2)(j + 2a + n - 2)!(j - n + 1)! 



2r(i + 2a + 3/2)(j-l)! 



(83) 



where ci 2 ™" 1 is an arbitrary constant. 



B. Semiclassical approximation to the C operator 

In this subsection we attempt a semiclassical calculation of the C operator. In such a calculation we expand Q as 
a series in powers of h: 



Q = Q +hQ 1 +h 2 Q 2 



(84) 



The semiclassical approximation terminates after the h term in this expansion. The ordering of powers of p and q in 
this expansion becomes unimportant because commuting p with q introduces additional powers of H. Furthermore, 
every factor of pq has dimensions of H, and thus only the first term in a sum needs to be kept. 

A recursive determination of the first-order solution Q\ in (84) arises as a natural simplification of the difficult 



problem that we formally solved in Subsec. (VA). To proceed, a dimensional analysis of the operators is required. 
Because of the commutator equation [q,p] — ih, we can assign the dimensions of q and p to be h 1 ^ 2 . With this 
convention the Hamiltonian for the shifted harmonic oscillator becomes 



H=kp 



1„2 , 12 



+ i 



h 1 ' 2 



For Qo in (191 to be the zeroth-order solution in h, we must introduce its explicit dependence on K; that is 

oo 

Qo — h ^ flfcJl-2fc,2fc- 



(85) 



(86) 



k=0 



Following the procedure illustrated in Sec. ( |IV[ ), we make the replacement Qo ~~ > A 4 Qo, where /i is a small dimensionless 
parameter. The operator Q\ in (74) admits dimensionless solutions in terms of the operators Rin only for a = in 



Introducing the explicit dependence on h, the operator Z 2 n becomes 



(|78|). In fact, for a = the series representation of both Z 2n in (|76|) and i?2„ in (|78|) can be drastically simplified. 

(87) 



7 _ fc » + i/ 2 y (-l)" +fc - 1 (fc-l)! T 

^2n — h 2_^t ~rZ. TTiTI -*-2fc,2fc-2n+l? 



k—n 



(n - l)!(fc - n)l 



while for the operators R 2n we get 



where the coefficients p^ n%> satisfy the recursion relation 



t> _ tn-1/2 „( 2 »)t 

tt-Zn — « Pfc J-2fc-l,2fc-2n+2) 



k=n— 1 



(2n) , 2fc + 1 (2n) _ 
Pfc+1 2(fc-n + 2) P,£ 2(n-l)!(*-n + 2)!' 



The general solution to (89) contains an arbitrary constant C n : 

(-l)»T(n + 1/2)[V5FW - r(fc - 1/2)] 



„(2") _ r 
Pk - °« 



/ 7 i : (n-l)!r(/c + l/2) 



(89) 



(90) 



With the choice C n = (— l) n 7r 1 / 2 r(n + l/2)/(n — 1)! the result in (90) is considerably simplified. The simplest 
first-order solution Q\ in the series (|84| for the operator Q is 



(-At 2 fi)" y> fc! 

(n-1)! ^ r(fc + l/2) 



T-2k-1.2k- 



2n+2- 



(91) 



n=0 v ' fe=n— 1 

This illustrates the nature of a semiclassical expansion for the operator Q. 
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VI. CONCLUSIONS 

The principal result in this paper is that while there is a unique bounded metric and C operator for the quantum 
harmonic oscillator, which is the simplest "PT-symmetric quantum theory, there is an infinite number of unbounded 
metric and C operators, and we have calculated them exactly. To produce these unbounded operators we have had 
to sum infinite series of singular operators (involving powers of 1/p) and have observed that the resulting sums are 
no longer singular. Of course, our summation procedure is at best only formal. However, we have verified our results 
by using dimensional summation procedures and have shown that the C operators that we have constructed satisfy 
exactly their defining equations. 

As anticipated in Ref. |31j . the properties of nonuniqucncss and unboundedness of the C operators are connected. 
There is a unique bounded C operator for the harmonic oscillator, namely V, and an infinite class of unbounded C 
operators. Interestingly, the unbounded metrics grow for large q like e q . This does not pose a serious problem if we 

2 

want to calculate matrix elements of eigenstates of Hq because in q space these states vanish like e~ q . Thus, any 
finite linear combination of eigenstates is an acceptable state in the Hilbcrt space associated with the unbounded 
metric. 

Finally, while we have performed formal summations of operators in this paper, we have justified our results by doing 
careful summation calculations that rely on analytic continuation. Our calculation are modeled on the dimensional 
continuation evaluations that are used to regulate divergent Feynman integrals. We conjecture that such techniques 
might be applied to generalize the notions of Cauchy sequences and completeness sums for Hilbert-space vectors. 
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